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Abstract 

In this paper we discuss massive gravity in Minkowski space via gravita- 
tional Higgs mechanism, which provides a non-perturbative definition thereof. 
Using this non-perturbative definition, we address the issue of unitarity by 
studying the full nonlinear Hamiltonian for the relevant metric degrees of 
freedom. While perturbatively unitarity is not evident, we argue that no 
negative norm state is present in the full nonlinear theory. 
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1 Introduction and Summary 



The gravitational Higgs mechanism gives a non-perturbative and fully covariant 
definition of massive gravity [Il[2l[3lll[5l[6l[71|8l|9l[l0l[IIl[l2l[l3l[Tll[T^ The 
graviton acquires mass via spontaneous breaking of the underlying general coor- 
dinate reparametrization invariance by scalar vacuum expectation values]^ In this 
paper, motivated by the results of [12] for the gravitational Higgs mechanism in de 
Sitter space, we study the issue of unitarity of the gravitational Higgs mechanism 
in the Minkowski background within the setup of [5]. 

Within the perturbative framework, unitarity requires that, in order not to prop- 
agate a negative norm state at the quadratic level in the action, the graviton mass 
term be of the Fierz-Pauli form [52] • Furthermore, higher order terms should be 
such that they do not introduce additional degrees of freedom that would destabi- 
lize the background, and difficulties in achieving this to all orders in perturbative 
expansion have been known for quite some time (see, e.g., jlQ], [27], and references 
therein). In this regard, in order to circumvent the aforesaid difficulties within the 
perturbative expansion, [SB] proposed an order-by-order construction, albeit in a 
special decoupling limit, such that higher-than-second-order time derivatives in the 
equations of motion are absent. 

Here we propose a different approach. Following fT2], our key observation is that 
perturbation theory appears to be inadequate, among other things, for the pur- 
poses of addressing the issue of unitaritylfl Thus, while the theories of [S] reproduce 
the Fierz-Pauli action at the quadratic level, according to [13], at higher orders in 
perturbative expansion they do not reduce in the decoupling limit to the theories 
studied in [36] (and, according to [H], the same holds for the models discussed in 
[9]). However, since the definition of [5] is intrinsically non-perturbative, we can test 
the stability of the Minkowski background in the full nonlinear theory. In this paper 
we perform a non-perturbative Hamiltonian analysis for the relevant metric degrees 
of freedom and argue that the full nonlinear theory appears to be free of ghosts. 

Our main result is that, in the gravitational Higgs mechanism, non-perturbatively, 
the Hamiltonian appears to be bounded from below in the Minkowski backgrounds 
We argue that this is indeed the cas^ by studying full nonlinear Hamiltonian for 
the relevant conformal and helicity-0 longitudinal modes with no spatial dependence, 

^For earlier and subsequent related works, see, e.g., [ll[T2l[ia[ia[20l[2Tl[2a[23[2l[25l[2i 
[27l[2i[2S[30l[ni[3l[33l32|3S|3i|32l[3^ references therein. 

^For instance, in the de Sitter case, this becomes evident from the fact that in the full nonlinear 
theory no enhanced local gauge symmetry (or ghost) is present for any value of the Hubble param- 
eter, while its appearance in the perturbative framework already at the quadratic order appears 
to be a mere artifact of linearization jl2| . 

^We arrived at the same conclusion in [12] in the de Sitter case. 

^Here we should note that in this regard we only analyze in detail the simplest example with 
higher (namely, four) derivative couplings in the scalar sector (^5)) : however, based on the fact that 
in the gravitational Higgs mechanism diffeomorphisms are broken spontaneously, we believe our 
conclusions should hold in the general case as well, albeit we do not have a proof of this statement. 
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which is the dimensionally reduced diagonal Ansatz of [12]. We show, however, that 
within the same Ansatz, depending on the choice of the field parametrization, a 
perturbative decomposition in terms of the conformal and helicity-0 modes invari- 
ably leads to a) equations of motion with higher-than-second-order time derivatives 
in agreement with the claim of |36] or b) a ghost already at the quadratic level. 
We emphasize that the results of [361 [IH EH] are obtained in the context of an 
intrinsically perturbative field parametrization, which parametrization does not ap- 
pear to posses a non-perturbative generalization, and all "no-go" results stemming 
therefrom appear to be mere artifacts of perturbative expansion. 

We also revisit the gravitational Higgs mechanism in the simplest case with no 
higher derivative couplings in the scalar sector, first discussed in [2], which does 
not correspond to the Fierz-Pauli mass term at the quadratic level. Nonetheless, 
non-perturbatively, even this case appears to "resum" into a theory with a posi- 
tive definite Hamiltonianj^ We argue that this is the case using our dimensionally 
reduced diagonal Ansatz. We also reproduce the same result using the full Hamil- 
tonian analysis of [TTj, thereby validating our Ansatz. 

The essence of our results is well-illustrated by the following simple "toy" exam- 
ple: 



where T-L is the Hamiltonian, p is the canonical conjugate momentum, and is a 
parameter. Non-perturbatively, this Hamiltonian is positive definite. On the other 
hand, in the "weak-field" , small p"^ approximation, perturbatively there is a ghost as 
the kinetic term proportional to has a wrong sign: Ti = fi—p"^ /2fj,+. . . ; simply put, 
in the above example the weak-field approximation is invalid in the regime where 
the fake perturbative "ghost", which is merely an artifact of linearization, would 
destabilize the background. The same appears to be the case in the gravitational 
Higgs mechanism examples we study in this paper. 

To summarize, non-perturbatively, the gravitational Higgs mechanism appears to 
be free of ghostsjj. 

The rest of the paper is organized as follows. In Sections 2 and 3 we discuss the 
gravitational Higgs mechanism in Minkowski background, which results in massive 
gravity with the Fierz-Pauli mass term for the appropriately tuned cosmological 
constant. In Section 4 we derive the Hamiltonian for the relevant metric modes 

®So does a continuous set smoothly connecting this case to the aforesaid Fierz-Pauh case. 

^However, in this paper we do not attempt to address the question of whether there is any 
superluminal propagation of signals or the related issue of causality. In this regard, we emphasize 
that the recent "no-go" results of [H] are obtained in the context of the aforesaid intrinsically 
perturbative field parametrization, directly rely on the results of [35], and do not apply to the full 
non-perturbative definition of the gravitational Higgs mechanism. To see if there is any superlu- 
minal propagation of signals in the full non-perturbative theory, it appears that one might have to 
develop some new non-perturbative methods, which is clearly beyond the scope of this paper. The 
non-Fierz-Pauli model of [2] is the "least non-perturbative" and might provide a fruitful testing 
ground in this context. 
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and show that it is bounded from below. In Section 5 we revisit the simplest case 
without higher derivative couplings in the scalar sector and show that there too the 
Hamiltonian is positive definite, in spite of the quadratic truncation not being of 
the Fierz-Pauli form. 



2 Minkowski Solutions 

Consider the induced metric for the scalar sector: 

Ymn = ZabVm<P^Vn(P'' . (2) 

Here M = 0, . . . , (D — 1) is a space-time index, and A = 0, . . . , {D — 1) is a global 
index. We will choose the scalar metric Zab to be the Minkowski metric: 

Zab = Vab ■ (3) 

Let 

Y ^ YmnG'''' . (4) 
The following action, albeit not the most genera will serve our purpose here: 

Sy = M^-^ j d^x^/^ [R - V{Y)] , (5) 

where a priori the "potential" V{Y) is a generic function of Y . 
The equations of motion read: 

v*^ (v'(r)VM0^) = , (6) 

Rmn - \GmnR = V'{Y)Ymn - \GmnV{Y) , (7) 

where prime denotes derivative w.r.t. Y . Multiplying ([6]) by Zab^ s'P^ and con- 
tracting indices, we can rewrite the scalar equations of motion as follows: 



d 



'M 



\/ZGV'iY)G''''YNs] - ^V^V'iY)G^'''dsYMN = . (8) 



""^"One can consider a more general setup where the scalar action is constructed not just from 
Y, but from Ymn, Gmn and eMo...Mr>_i, see, e.g., [51 |H1 HI UHl HI]- However, a simple action 
containing a scalar function V{Y) suffices to capture all qualitative features of gravitational Higgs 
mechanism. In particular, if this function is quadratic as in (j23l) . the cosmological constant A must 
be negative in the context of the Minkowski background (but not in the de Sitter case - see [H]); 
however, generically there is no restriction on A, which can be positive, negative or zero even in the 
context of the Minkowski background, once we allow cubic and/or higher order terms in V{Y), or 
consider non-polynomial V{Y). As a side remark, let us note that no choice of polynomial V{Y) 
constructed from Y only reduces in the decoupling limit to the theories studied in [3 6) (c/. |14j). 
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Since the theory possesses full diffeomorphism symmetry, ([8]) and ([7]) are not all 
independent but linearly related due to Bianchi identities. Thus, multiplying ([7]) by 
y/—G, differentiating w.r.t. and contracting indices we arrive at (IE]). 
We are interested in finding solutions of the form: 

(P^ = m 6^M , (9) 
Gmn = Vmn , (10) 

where m is a mass-scale parameter. The equations of motion ([7]) imply that 

Y, = Dm^ (11) 

is the solution of the following equation 

V{Y,) = I Y,V'{Y,) , (12) 
which determines the mass scale m. 



3 Massive Gravity 

In this section, following [5], we study linearized fluctuations in the background 
given by ([9]) and (ITOl) . Since diffeomorphisms are broken spontaneously, the equa- 
tions of motion are invariant under the full diffeomorphism invariance. The scalar 
fluctuations (p^ can therefore be gauged away using the diffeomorphisms: 

= VM<p^e' = m . (13) 

However, once we gauge away the scalars, diffeomorphisms can no longer be used to 
gauge away any of the graviton components Hmn defined as: 

Gmn = Vmn + hjifw ■ (14) 

Moreover, we will use the notation h = rj^'^'^hMN- 
After setting (p^ = 0, we have 

Ymn = m^VMN , (15) 
Y = YmnG^^^ = m^[D-h + ...] = Y,-m^h + ... , (16) 

where the ellipses stand for higher order terms in Hmn- 

Due to diffeomorphism invariance, the scalar equations of motion ([6]) are related 
to ([7]) via Bianchi identities. We will therefore focus on ([7]). Let us first rewrite it 
as follows: 

Rmn — -GmnR = 

[VmnV'{Y) - GmnV'{Y,)] - ^Gmn [V{Y) - V{Y,)] . (17) 
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Linearizing the r.h.s. of this equation, we obtain: 

1 

Rmn — -^GmnR = [VMNh — C^MAf] + • • • , (18) 

where 

M2 = mV(K) - 2m^V"(K) , (19) 
(M^ = 2mV(y;) . (20) 

This corresponds to adding a graviton mass term of the form 

- ^ [ChMNh'''' - h'] (21) 

to the Einstein-Hilbert action, and the Fierz-Pauh combination corresponds to tak- 
ing C = 1. This occurs for a special class of potentials with 

V'{Y,) = -|i;V"(K) . (22) 

Thus, as we see, we can obtain the Fierz-Pauli combination of the mass term for 
the graviton if we tune one combination of couplings. In fact, this tuning is nothing 
but the tuning of the cosmological constant - indeed, (122|) relates the cosmological 
constant to higher derivative couplings. 
Thus, consider a simple example: 

V = A + Y + XY^ . (23) 

The first term is the cosmological constant, the second term is the kinetic term for 
the scalars (which can always be normalized such that the corresponding coefficient is 
1 by normalizing the scalars (p^ accordingly) , and the third term is a four-derivative 
term. We then have: 

Y* = , ^ , A-^ , (24) 

2(D + 2) ' ^ ^ 

which relates the mass parameter m to the higher derivative coupling A: 

= YJD = A~^ , (25) 

and the graviton mass is given by: 

= A-i . (26) 

(D + 2)2 ^ ' 

Note that we must have A < 0. Moreover, we have: 

_ I^^ + 4D-8 

4(D + 2)2 ^ • ^27) 

So, the cosmological constant in this case must be negative, which is due to the 
choice of the potential ( 12^ : however, as we already emphasized above, for generic 
choices of the potential there is no restriction on the cosmological constant, which 
can be positive, negative or zero. 
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4 Is there a Ghost? 



The purpose of this section is to argue that the full nonlinear theory of massive 
gravity in Minkowski space via gravitational Higgs mechanism is free of ghosts. We 
will do this by studying the full nonlinear action for the relevant modes, which we 
identify next. In particular, we will argue that no negative norm state is present for 
these modes. 

Let us note that, once we gauge away the scalars, the full nonlinear action 
becomes: 

~ ~ ' ^ ~ ' (28) 



where V{C) = V{w?Q). 

To identify the relevant modes in the full nonlinear theory, let us note that in the 
linearized theory the potentially "troublesome" mode is the longitudinal helicity-0 
mode p. However, we must also include the conformal mode u as there is kinetic 
mixing between p and u. In fact, p and u are not independent but are related via 
Bianchi identities. Therefore, in the linearized language one must look at the modes 
of the form 

huN = Vmn + Va/ VatP . (29) 
Furthermore, based on symmetry considerations, namely, the SO{D — 1) invariance 
in the spatial directions]^ we can focus on field configurations independent of spa- 
tial coordinates [12]. Indeed, for our purposes here we can compactify the spatial 
coordinates on a torus T^~^ and disregard the Kaluza-Klein modes. This way we 
reduce the D-dimensional theory to a classical mechanical system, which suffices for 
our purposes here. Indeed, with proper care (see [I2j), if there is a negative norm 
state in the uncompactified theory, it will be visible in its compactified version, and 
vice- versa. 

Let us therefore consider field configurations of the form: 

G*^^ = diag(^?(t)r/°°,/(t)r/^^) , (30) 

where g{t) and f{t) are functions of time t only. The action ( 128|) then reduces as 
follows: 



SG = -t^ j dt g-^^f-"^ {jgU' + V{g + Q)] 



(31) 



where 



K = M^-^Wd^i , (32) 

^^{D-l){D-2) , (33) 

U^^dMf), (34) 

n^iD-l)f , (35) 



^^Indeed, negative norm states cannot arise from purely space-like components or spatial deriva- 
tives, and are due to time-like components and/or time derivatives. 
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and Wd-i is the volume in the spatial dimensions {i.e., the volume of T^~^). Note 
that is a Lagrange multiplier. The goal is to integrate out g and obtain the 
corresponding action for /. It is then this action that we should test for the presence 
of a negative norm state. 

The equation of motion for g reads: 



V{g + fi) - 2gV\g + fi) = -fgU^ 



(36) 



The following discussion can be straightforwardly generalized to general V . However, 
for our purposes here it will suffice to consider quadratic V corresponding to 
We then have: 



1 + 2\m^n + ^U^ 



A 



+ [l + Xrn^n] > = . 



(37) 



We can therefore express g in terms of / and (9tln(/): 
6Xm^g = 



1 + 2Am2fi + At f/^ 



1 + 2Am2^] + ^ f/2 



+ I2\m? + + Xmm] \ , (38) 



where the branch is fixed by the requirement that g = 1 when / = 1. Substituting 
the so expressed g into ( l3T]) . we obtain an action which is a nonlinear functional of 
/ and 9i ln(/). 

For our purposes here it is more convenient to work with the canonical variable 
g, where 



q = ln(/) , 



and the action reads: 



G 



dt L = —K I dt g 2e 



^gU^ + V{g + Q) 



(39) 
(40) 

(41) 



(42) 



where L is the Lagrangian. This action corresponds to a classical mechanical system 
with a lagrange multiplier g. 

Next, the conjugate momentum is given by 



P 



dL 
d{dtq) 



X 



X <i ^g'HlU^ - ^9'Hv{g + Q)+ g-Hv'{g + Q) + g'^^U 



(43) 



where 



9 = 



dg 
d{drq) 



Using (136!) . (143|) simplifies to 

p = —ne 2 
and the Hamiltonian is given by 



U=pdtq-L = -ng-^e — ^'^ -fgU^ - V{g + fi) 

We can now see if this Hamiltonian is bounded from below. 
First, using f p6|) . we have: 



n = 2Kgh~^'^V\g + fi) = 2m^Kgh-^'^ [l + 2\m^{g + Vt 

Using (l38l) . we can rewrite this Hamiltonian as follows: 

2 



(44) 



(45) 



(46) 



1 D-l . 



U = -m^ng^e-^'^ [X - Z] , 
3 



(47) 



(4^ 



where 



X = 



1 + 2\m?VL + ^ f/2 



Z = ^U^ - AXm^^l - 2 . 



+ 12Am2 <; A + [1 + \mm I , (49) 

(50) 



The presence of a ghost would imply that the Hamiltonian is unbounded from below 
for large values of f/^ (recall that U"^ contains the "kinetic" term). However, it is 
not difficult to show that this Hamiltonian suffers from no such pathology. Indeed, 
we can rewrite it as follows: 



2 2 1 -D-i „ X'^ — Z'^ 

3 ^ X + Z 



2rn?Kg2e 2 1 
which in the large f/^ limit reads: 



1 4AA - 1 + 2(7/m2) U'^ [1 + 2\m'^n] 



X + Z 



n = 2m\gh-^'' [l + 2\m^n] + 0{1/U^ 
Furthermore, from fl38ll we have 



X + Q 



X + Q 



(51) 



(52) 



(53) 



where 
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Q^^U' + 2Xm'n + 1 . (54) 



So, in the large f/^ hmit the Hamiltonian actually vanishes. 

Note that the above argument implicitly assumes that Q is bounded from above. 
This is indeed the case as g must be at least non-negative, which implies that 

AA + Xm^Q [1 + Xm^Q] < , (55) 

and Q is bounded as follows (note that we must have > 0): 

and we must further have 

A>^. (5T) 

which is always satisfied due to (127|) as A < 0. 

Thus, as we see, there appears to be no ghost in the full nonlinear theory. The key 
ingredient here is the parametrization of the conformal and helicity-0 longitudinal 
modes. We have been working with (1301) . while the results of [Ml 111] apply to (1291) . 
The difference between the two is that fl30l) has no derivatives. Thus, at the quadratic 
order, the second derivatives introduced by the parametrization corresponding to 
f l29|) can be integrated by parts to arrive at an action containing only first derivatives 
of uj and p. However, we have explicitly checked that already at the cubic level the 
second derivatives introduced by the parametrization corresponding to ( 129|) cannot 
be integrated by parts, so the resulting action invariably includes terms with second 
derivatives of p. This is clearly problematic already at the cubic level and suggests 
that the parametrization corresponding to ( |29|) cannot be used beyond the linearized 
approximation. Indeed, a nonlinear completion of (!29|) is given by: 

where the covariant derivative is defined w.r.t. the metric timn (this choice does 
not affect our discussion here). Note, however, that such a parametrization of the 
metric is rather problematic in the context of the full nonlinear theory fl28|) as it 
introduces higher derivative terms in u, which should therefore not be used as the 
canonical variable in the full nonlinear theory. This suggests that our parametriza- 
tion ( 130|) is indeed adequate (while our results here suggest that the perturbative 
field parametrization underlying the results of [H] does not appear to be adequate 
in the full non-perturbative theory). 

5 The Linear Potential Revisited 

In this section we consider the linear potential 

V{Y) = K + Y . (59) 
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This example was discussed in [2j and worked out in detail in the Hamiltonian for- 
malism in [TT]. The potential (15^ does not satisfy the condition (122]) and, therefore, 
does not lead to the Fierz-Pauli action at the quadratic order; consequently, it propa- 
gates a negative norm excitation at that order in perturbation theory. Nonetheless, 
we will show that the full non-perturbative Hamiltonian is bounded from below. 
This is shown in two ways: In the dimensionally reduced parametrization of fl30l) . 
and in the full Hamiltonian analysis of [Hj. This supports our conclusion that the 
Ansatz (1501) is fully adequate for analyzing unitarity. 

Recalling the definitions of Section HJ we have V{x) = A + rn^x, and the g 
equation of motion gives 

g = {A + m^Q)/{m^ + ^U^), (60) 

thus, the Hamiltonian reads: 

n = 2Km^ e^^^^i/' , (61) 

which is in general positive definite (assuming g > 0, which is the case for / > 
-A/{D-l)m'^ = {D-2)/{D-l), where we have used and (^). Therefore, as 
for the quadratic potential of the previous section, it appears that non-perturbatively 
there is no ghost. Furthermore, let us point out that perturbatively there is a ghost, 
since in a "weak-field" , small f/^ expansion the kinetic term proportional to has a 
wrong sign. This is similar to instructive toy examples discussed in Section 1 hereof 
as well as Appendix A of [T2|. 

We arrive at the same result in full generality by introducing the lapse constraint 
into the gauge-fixed Hamiltonian derived for this potential in [TT] . The result reads, 
in their notation, 

which is positive definite, and reduces to ( ET\i for the Ansatz (15U]) since, in that case, 
nf^ = 0, = g and det h = l/f^-\ 

In fact, there is nothing "special" about the linear potential fl59|) in terms of 
unitarity. The discussion in Section 4 is completely independent of fl22|) (with the 
exception of the reference to ( 1271) immediately following ( 1571) . which, however, does 
not affect our discussion here - see below), which is the condition corresponding to 
having the Fierz-Pauli term at the quadratic order. In particular, the discussion in 
Section 4 and its conclusion that non-perturbatively there appears to be no nega- 
tive norm state are valid for all values of A, including those smoothly interpolating 
between (p7|) (which corresponds to the Fierz-Pauli term) and A = 0, which corre- 
sponds to the linear potential (15^ as long as (ITTI) is satisfied along the interpolation 
path (and this condition does not pose an obstruction to such interpolation). Simply 
put, non-perturbatively there appears to be no negative norm state for a continuous 
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family of models arising out of the gravitational Higgs mechanism, and the linear 
potential point and the Fierz-Pauli point are not any special in this regard. 

The fact that non-Fierz-Pauli points appear to be unitary deserves further elab- 
oration. Thus, in the gravitational Higgs mechanism we start with massless gravity 
with D[D — 3)/2 degrees of freedom coupled to D scalars, one of which is time-like. 
As was explained in [5], at the Fierz-Pauli point fl22|) the time-like scalar does not 
propagate, so in the scalar sector we have only {D — 1) propagating space-like de- 
grees of freedom, which are eaten in the gravitational Higgs mechanism producing 
{D + !)(£'- 2)/2(= D{D - 3)/2 + {D - 1)) degrees of freedom for the massive 
graviton, all of which are unitary. What about the non-Fierz-Pauli points? Here we 
have an extra time-like scalar degree of freedom, which in the perturbative language 
in the Higgs phase results in a ghost-like propagating degree of freedom, namely, 
the trace of the graviton h = m- However, non-perturbatively the Hamiltonian is 
positive-definite and there is no ghost, while if the Hamiltonian is expanded pertur- 
batively, as we saw above, a ghost invariably appears. Simply put, the weak-field 
approximation is invalid in the regime where the fake perturbative "ghost", which 
is merely an artifact of linearization, would destabilize the background!^ 
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